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Abstract 

In the light of the Sturm-Liouvihe theorem, the Levinson theorem for the 
Schrodinger equation with both local and non-local cylindrically symmetric po- 
tentials is studied. It is proved that the two-dimensional Levinson theorem holds 
for the case with both local and non-local cylindrically symmetric cutoff po- 
tentials, which is not necessarily separable. In addition, the problems related 
to the positive-energy bound states and the physically redundant state are also 
discussed in this paper. 
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1. Introduction 



The Levinson theorem [1], an important theorem in the scattering theory, estab- 
hshcd the relation between the total number of bound states and the phase shifts at 
zero momentum. The Levinson theorem has been proved by several authors with dif- 
ferent methods, and generalized to different fields [2-9] . Generally speaking, three main 
methods are used to prove the Levinson theorem. One [1] is based on the elaborate 
analysis of the Jost function. This method requires nice behavior of the potential. 
The second, based on the Green function method [5], expounds that the total number 
of the physical states, which is infinite, is proved to be independent of the potential 
and the number of the bound states is the difference between the infinite numbers of 
the scattering states without and with the potential. The third method proves the 
Levinson theorem by the Sturm-Liouville theorem [6-8] . This simple, intuitive method 
is easy to be generalized. Some obstacles and ambiguities, which may occur in other 
two methods, disappear in the third method. We have succeeded in dealing with the 
nonrelativistic and relativistic problems in two dimensions by this way [10-11]. 

The reasons why we write this paper are that, on the one hand the Levinson theorem 
in two dimensions has been studied in experiment [12] as well as in theory [13,10,11] be- 
cause of the wide interest in lower-dimensional field theories and other modern physics 
[14-20], on the other hand the Levinson theorem for non-local interactions in two di- 
mensions has never appeared in the literature. 

This paper is organized as follows. In Sec. 2, we estabhsh the Sturm-Liouville 
theorem for non-local interactions in two dimensions. The Levinson theorem for this 
case will be set up in Sec. 3. Some problems related to the positive-energy bound 
states and the physically redundant state will be studied in Sec. 4 and 5. 

2. The Sturm-Liouville Theorem 

Throughout this paper h = 1 and the mass = 1/2 are employed for simplicity. 
Consider the Schrodingcr equation with a local potential V{r) and a non-local potential 
U{r,r'), where both potentials are cylindrically symmetric 

fid d 1 (9M 
Hij{r,^) = - [-—r— + -^]^{r,^) + V{r)4j{r,^) 



+ 



J U{r, r')S{^' - ^)ip{r', ip')r' dr' dip 



= Elixir, V)- 

The mesonic theory of nuclear forces points out the interaction between two nucle- 
ons is local at great distances but becomes non-local if the two nucleons come close. To 
simplify the expression, we assume, following Martin [4], that the non-local potential 
U (r, r') is real, symmetric, vanishing at large distance, and not too singular at the 
origin [21], 

C/(r, r') = C/(r', r) 

yU{r,r')\ ~ 0, at r ~ 0, (2) 
C/(r, r') = 0, when r > Tq. 
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As usual, the local potential V{r) should not be too singular at the origin and at 
infinity. We assume that V{r) satisfies 

|rV(r)| - 0, at r ~ 0, , . 

V{r) = 0, when r > tq. 

The first condition is necessary for the nice behavior of the wave function at the origin 
[1], and the potential with the second condition is called a cutoff one, namely, it is 
vanishing beyond a sufficiently large radius Tq. It was proved that the tail of the local 
potential at infinity will not change the essence of the proof [10] if it decays faster than 
at infinity. Under this assumption, the integral range in (1) is, in fact, from to 
To, and the equation in the region [tq, oo) becomes that for the free particle. 
Introduce a parameter A for the potentials 

V{r, A) = AT/ (r) , U(r, r'. A) = XU(r, r') . (4) 

As A increases from zero to one, the potentials V{r, A) and U (r, r'. A) change from zero 
to the given potentials V{r) and U{r,r'), respectively. 
Owing to the symmetry of the potentials, letting 

V^(r, 0, A) - r-V'i?E^(r, A)e±^"^^ m = 0, 1, 2, . . . , (5) 

we obtain the radial equation 



— i?£;^(r, X)+{E- V{r, A) j i?s^(r. A) 



(6) 



where ±m and E denote the angular momentum and energy, respectively. Since the 
radial function REm{f-i A) is independent of the sign of the angular momentum, we only 
discuss the case with non- negative m. 

We are going to solve the radial equation (6) in two regions [0,ro] and [ro, oo), and 
match two solutions at Tq. Since the Schrodinger equation is linear, the wavefunction 
can be multiplied by a constant factor. Removing the effect of the factor, we only need 
one matching condition at Tq for the logarithmic derivative of the radial function: 

\RErn{r,\) dr J [i?E„,(r,A) dr J ^^^^^ 

We now turn to the Sturm-Liouville theorem. Denote by REmir, A) a solution of 
(6) with the energy E. Multiplying the equations for REmit^,^) and REm{f-i^) by 
REmir, A) and REmij, A), respectively, and calculating their difference, we have 

^ \REm{r, A) — REm{r, A) — j +{E- E)REm{r, X)REm{r, A) 



(8) 



^/TREUr, A) J Uir,r',X)RErnir',\)V?dr' 
-V^Em{r,X) j U{r,r',X)REm{r'A)V?dr'. 
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According to the boundary condition, both REmif^, A) and REm{r, A) go to zero when 
r tends to zero. Integrating (8) over the variable r in the region [0, Tq] and noting the 
symmetric property of U (r, r') , we have 

1 ^R,^^r,xf^^!^-R,^(r,xf''-'"^^^''^'^ 



E — E [ ' dr ' dr 

fro 



pro 

- / REm{r', X)REm{r\ X)dr'. 
Jo 



'0 

Taking the hmit, we obtain 

dA^{E,X) ^ d ( 1 dREm{r,X) \ 
dE - dE\RErn{rA) dr J 

= - REm{ro, A)-2 / REm{r', Xfdr' < 0. 

JO 

Similarly, from the boundary condition that when E <Q the radial function REm{f-i A) 
tends to zero at infinity, we have 



d ( 1 dREm{r,\) \ 

dE \REm{r,X) dr J_^^ " L 



/•OO 

REmiro, A)-2 / REm{r', xfdr' > 0. (10) 



Therefore, when E both sides of the matching condition (7) are monotonic with 
respect to the energy E. As energy increases, the logarithmic derivative of the radial 
function at Tq— decreases monotonically, but that at ro+ increases monotonically. This 
is an expression for the Sturm- Liouville theorem [22] . 

3. The Levinson Theorem 

The establishment of the Levinson theorem for the case with both local and non- 
local cylindrically symmetric potentials is similar to that for the case with only a local 
potential. 

In solving the radial equation (6) in the region [0, Tq], only one solution is convergent 
at the origin. Thus, for the given potentials the logarithmic derivative A^{E, X) is 
determined in principle. For example, for free particle (A = 0) we have 



]^—^Jjn{kr), when E>0 and k — Ve ^ ^ 

when E <0 and k = y/^, 



REm{r, 0) 



where the factor in front of the radial function REm{i^) is not important. The solution 
REm^T-i 0) given in (11) is a real function. The logarithmic derivative at tq— for £^ < 
is 



dREm{r,0) 



REmir, 0) dr 



r=ro- 



[ kJ^ikro) 1 ^^^^ ^ ^ 
Jm{kro) 2ro 

VlM!^_l_ when E<0 

(12) 
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In the region [vq, oo), we have V{r) = U{r,r') = 0. For £" > 0, there are two oscil- 
latory solutions to (6). Their combination can always satisfy the matching condition 
(7), so that there is a continuous spectrum for £^ > 0. 



REm{r,X) = J— -{cos r]m(k,X)Jm{kr) -sin r]m{k,X)N^(kr)} 

' / X (13) 

~ cos \ kr — + 7]m{k, A) j , when r — > oo, 

where N^{kr) is the Neumann function. The phase shift rj^ik. A) is determined by the 
matching condition (7) 

. . _ Mkro) A,„(£. A) -/,■./;,(/,■/■„)/./,„ (/,ro)-l/2ro 
iann^[K,A) ^^^^^^^ ■ j^^^E,\) - kN'^{kr,)/N^{kr,) - l/2ro' ^ ^ 

rim{k) =ri^{k,l), (15) 

where the prime denotes the derivative of the Bessel function, the Neumann function, 
and later the Hankel function with respect to their argument. Although the radial 
equation (6) in the region [tq, oo) is independent of A, the solution REm{f-i A) and the 
phase shift i]„i{k, A) do depend on A through the matching condition (7). 

The phase shift r]m{k, A) is determined from (14) up to a multiple of vr due to the 
period of the tangent function. Levinson determined the phase shift rim{k) with respect 
to the phase shift 77^(00) at the infinite momentum. For any finite potential, the phase 
shift ?7m(oo) will not change and is always equal zero. Therefore, Levinson's definition 
for the phase shift is equivalent to the convention that the phase shift ri.,n{k,\) is 
determined with respect to the phase shift 'qm{k, 0) for the free particle, where rim{k, 0) 
is defined to be zero 

r]m{k, 0) = 0, where A = 0. (16) 

There is some ambiguity for rim{oo) when a bound state with a positive energy occurs 
(see Sec. 4). However, as far as the Levinson theorem is concerned, the latter con- 
vention is more convenient. We prefer to use this convention where the phase shift 
f]m{k) is determined completely as A increases from zero to one. It is the reason why 
we introduce the parameter A. 

For E <0 there is only one convergent solution at infinity 



i?Em(r) = e^("^+^)'^/^/^i^«(^«;r)~e-'^^ when r ^ 00. (17) 



where H^^^x) is the Hankel function of the first kind. Thus, the matching condition 
(7) is not always satisfied. When the matching condition (7) is satisfied, a bound state 
appears at this energy. It means that there is a discrete spectrum for £■ < 0. 
From (17) we have 

r 1 dREm{r,0)\ iKH^^^iKroY 1 



REmir,0) dr j^^^^^ Hin\inro) 2ro 

f (-m + l/2)/ro = Pm when E — >0 
1 —K, ~ —00 when E — >■ —00. 

(18) 
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On the other hand, if V{r) = U{r,r') = 0, we obtain from (12) 



[ 1 dREm{r,0) 
\REm{r,0) dr 
iKJ!^{iK,ro) 1 _J (m + l/2)/ro when E — >0 



Jmiii^ro) 2ro 1 fi; oo when E — > -oo. 



It is easy to see from (18) and (19) that as energy increases from — oo to 0, there is 
no overlap between two variant ranges of the logarithmic derivatives at two sides of 
ro such that there is no bound state when A = except for S wave where there is a 
half bound state at ii^ = 0. The half bound state will be discussed at the end of this 
section. 

If Am{0, A) decreases across the value pm = (— ^+l/2)/ro as A increases, an overlap 
between two variant ranges of the logarithmic derivatives at two sides of Tq appears. 
Since the logarithmic derivatives of the radial function at Tq— decreases monotonically 
as the energy increases, and that at ro+ increases monotonically, the overlap means that 
there must be one and only one energy where the matching condition (7) is satisfied, 
i.e. a scattering state changes to a bound state. 

As A increases, a zero point in the zero energy solution RQm{r, A) may comes through 
ro- In this process ^^(0, A) may decrease to the negative infinity, jumps to the positive 
infinity, and decreases again, or vice versa. It is not a singularity. If A^(0, A) decreases, 
through the jump at infinity, again across the value p^, another bound state appears. 

As A increases from zero to one, each time A^{0, A) decreases across the value 
p„i, a new overlap between the variant ranges of two logarithmic derivatives appears 
such that a scattering state changes to a bound state. Conversely, each time A,n{0, A) 
increases across the value pm, an overlap between those two variant ranges disappears 
such that a bound state changes back to a scattering state. The number of bound 
states rim is equal to the times that ^^(0) decreases across the value pm as A change 
from zero to one, subtracted by the times that ^^(0) increases across the value p^. In 
the following, we will show from (14) that this number is nothing but the phase shift 
rim{0) at zero momentum divided by tt. 

It is easy to see from (14) that the phase shift r]m{k, A) increases monotonically as 
the logarithmic derivative Ajn{E) decreases 



dr]mik,X) 



dAm{E,X) 



— 2 ^ 0- (20) 



TT {2roAm{E)Nm{kro) - 2kroN^{kro) - Nm{kro)y 



The phase shift ?7m(0. A) is the limit of the phase shift rim{k,\) as k tends to zero. 
Therefore, what we are interested in is the phase shift rim{k,X) at sufficiently small 
momentum k, A; -C I/tq. For the small momentum k we obtain from (14) 
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tanr]m{k, A) ~ 

-7r(A;ro)2™ A^(0, A) - (m + l/2)/ro 



when m > 2 



m — l)(2m — l)y 

7r(A:ro)^ yl^(0, A) - 3/(2ro) 

4 ■ ^„(0,A)-c2^2_^^(i + 2(fcro)2 1og(fcro)) w en m 

TT ^^(0,A)-c2A;2-po(l-(A;ro)') 



when m = 0. 



(21) 

In addition to the leading terms, we include in (21) some next leading terms, which is 
useful only for the critical case where the leading terms are canceled with each other. 

First of all, it can be seen from (21) that tan?7m(/^, A) tends to zero as k goes to 
zero, i.e. r/m(0,A) is always equal to the multiple of vr. In other words, if the phase 
shift ri„,.{k,\) for a sufficiently small k is expressed as a positive or negative acute 
angle plus nir, where n is an integer, its limit rim{0, A) is equal to nvr. It means that 
77^(0, A) changes discontinuously. By the way, in three dimensions, the phase shift at 
zero momentum of S wave may have an additional 7r/2 when the half bound state 
occurs. 

Secondly, since the phase shift ri„,,{k,X) increases monotonically as the logarithmic 
derivative ^^(-E', A) decreases, the phase shift at zero momentum 77^(0, A) will jump 
by TT if tan?7m(fc. A) at the sufficiently small k changes sign from positive to negative 
as Am{E, A) decreases, and vice versa. When A changes from zero to one continuously, 
each time Am{0, A) decreases from near and larger than the value pm to smaller than 
that value, the denominator in (21) changes sign from positive to negative and the 
remaining factor keeps positive, such that the phase shift at zero momentum rim{0, A) 
jumps by tt. Conversely, each time ^4^(0, A) increases across the value p^, the phase 
shift at zero momentum /^^(O, A) jumps by — tt. Therefore, the phase shift r)m{0)/'!r 
is just equal to the times A„i{0,X) decreases across the value pm as A increases from 
zero to one, subtracted by the times Am{0, A) increases across that value. Therefore, 
we proved the Levinson theorem for the Schrodinger equation in two dimensions for 
non-critical cases: 

?7m(0) = rimn. (22a) 

Thirdly, we should pay some attention to the case of m = 0. When Ao{0,X) 
decreases across the value po — l/(2ro), both the numerator and denominator in (21) 
change signs, but not spontaneously because the next leading terms in the numerator 
and denominator in (21) are different. It is easy to see that the numerator changes 
sign first, and then the denominator changes sign, i.e. ta.nrio{k, X) at small k changes 
firstly from negative to positive, then to negative again such that 770(0, A) jumps by tt. 
Similarly, when Am{0, A) increases across the value po, rjo{0, A) jumps by —tt. 

When A = and m = 0, the numerator in (21) is equal to zero, the denominator 
is positive, and the phase shift i]o{0) is defined to be zero. If Aq{E) decreases when A 
increases from zero, the numerator becomes negative firstly, and then the denominator 
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changes from positive to negative such that the phase shift ?7o(0) jumps by n and 
simultaneously a new bound state appears. 

Finally, we turn to discuss the critical cases where a half bound state occurs. If the 
logarithmic derivative A^(0, 1) is equal to the value pm, the following solution with 
zero energy in the region [ro, oo) will match this Am{0, 1) at ro 

Rom{r, 1) = r-+V2. 

It is a bound state when m > 2, but called a half bound state when m — 1, and 0. A 
half bound state is a zero energy solution of the Schrodinger equation which is finite 

but does not decay fast enough at infinity to be square integrable. We are going to 
discuss the critical case where ^^(0, A) decreases (or increases) and reaches, but not 
across, the value pm as A increases from a value, a little smaller than one, to one. For 
definiteness, we discuss the case where ^^(0, A) decreases and reaches the value pm as 
A increases to one. In this case a new bound state with zero energy appears for m > 2, 
but does not appear for m = 1 and 0. We need to check whether or not the phase shift 
r7m(0) increases an additional vr. 

It is evident to find that the denominator in (21) for m, > 2 has changed the sign 
from positive to negative as ^^(0) decreases and reaches the value pm, i.e. the phase 
shift r)m{0) jumps by tt and simultaneously a new bound state of zero energy appears. 

For m = the next leading term of the denominator in (21) is positive and larger 
than the term — c^A:^, such that the denominator does not change sign, i.e. the phase 
shift ?7m(0) docs not jump. It meets the fact that no new bound state appears. 

For 171 = 1 the next leading term of the denominator in (21) is negative such that 
the denominator does change sign, i.e. the phase shift rj^iO) jumps by tt as A^(0) 
decreases and reaches the value pi. However, in this case no new bound state appears 
simultaneously. 

The discussion for the cases where ^^(0) increases and reaches the value Pm is 
similar. Therefore, Levinson's theorem (22a) holds for the critical cases except for 
m = 1. In the latter case, Levinson's theorem for the Schrodinger equation with both 
local and non-local interactions in two dimensions becomes: 

^m(O) = {rim + 1) TT, whcu m = 1, and a half bound state occurs. (22b) 

As discussed above, it is found that the Levinson theorem holds without any mod- 
ification for the case where a non-local potential is included. 

4. Positive-energy bound states 

It is well known that, in the case with only a local interaction, the wavefunction and 
its first derivative would never vanish at the same point except for the origin, so there 
is no positive-energy bound state. However, in the case with a non-local interaction, 
Martin showed that the solution with an asymptotic form is not unique when the 
potential satisfies some conditions [4], i.e. there exists the positive-energy bound state 
with a vanishing asymptotic form. If a small perturbative potential is added such that 
the non-local potential satisfies the conditions, the positive-energy bound state will 
appear and the phase shift at this energy increases rapidly by almost n. This can be 
seen explicitly in the examples given by Martin [4] and Kermode [23]. 
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It was pointed out by Kermode that the inverse tangent function is not single- 
valued and it is physically more satisfactory to include a jump by tt to the phase shift 
at the energy Eq, where a positive-energy bound state occurs. Martin and Chadan 
[4,21] defined the phase shift to be continuous even at Eq so that an additional tt will 
be included into 5(0) — 5(oo) for each positive-energy bound state. This is their reason 
to modify Levinson's theorem by the term an where a denotes the number of positive- 
energy bound state. But according to the viewpoint of Kermode, no modification to 
the Levinson theorem is required. 

However, the phase shift at zero energy in our convention doesn't change, no matter 
which viewpoint is used, i.e. no matter whether the phase shift jumps or not at the 
energy with a positive-energy bound state. Therefore, the Levinson theorem (22) holds 
for the cases where positive-energy bound states may occur. 



5. Redundant state 



The resonating group model of the scattering of nuclei, or other composite sys- 
tems, derives an effective two-body interaction in which a non-local potential appears. 
There are some physically redundant states which describe Pauli-forbidden states for 
the compound system, and the physical two-body states must be orthogonal to these 
redundant states [24]. In the case of three dimensions, Saito [25], Okai, et al [26] and 
Englefield-Shoukry [27] proposed a simple non-local term which guarantees that the re- 
quired orthogonality, and verified that it was a good representation of the interactions. 
If there is just one redundant state represented by the real normalized wavefunction 
U{r), then the two-dimensional Saito's equation is 

^,RE^{r) +(e- V{r) - ^^^i^) Re^It) 



roo 

Jo 



U^(s)ds = 1, 



and 



roo 

E / Uir)RErnir)dr = 0. (24) 
Jo 

The solution of (23) satisfies the orthogonality constraint except for that of zero energy. 
Saito's non-local potential is separable. 

If the Schrodinger equation with only a local potential V{r) has a bound state with 
a negative —£ < 0, the corresponding wavefunction is denoted by ip{r): 



^^(r) - [V(r) + ^(r) = E^(r), ^^^^ 



poo 

Jo 



1. 



It is obvious that U{r) — -^{r) satisfies (23) with zero energy. Therefore, it is 
the so-called physically redundant state. As far as a mathematical equation (23) is 
concerned, the redundant state is one of the bound states with zero energy. 
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